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I/-) , Models of random walks are considered in which walkers are born at one location 
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and die at all other locations with uniform death rate. Steady-state distributions of 
random walkers exhibit dimensionally dependent critical behavior as a function of 
the birth rate. Exact analytical results for a hyperspherical lattice yield a second- 



order phase transition with a nontrivial critical exponent for all positive dimensions 
/ 2, 4. Numerical studies of hypercubic and fractal lattices indicate that these 
exact results are universal. Implications for the adsorption transition of polymers at 
curved interfaces are discussed. 
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In this letter we study random walk models in which walkers are born at one location and 
die everywhere else. We believe that our analytical results for the critical behavior exhibited 
by spherically symmetric random walks are universal and therefore provide insights into the 
dynamics governing dissipative systems that interact at an interface As an application, 
we discuss the adsorption transition of a polymer near an attractive boundary 0. 

Our spherically symmetric random walk model uses a lattice that drastically simplifies 
the treatment of a dissipative system interacting with a curved boundary. For example, we 
can model polymer growth near a hyperspherical boundary in any dimension D > by a 
random walk on a one-dimensional lattice. Such a lattice faithfully represents the spatial 
entropy of the system. While this lattice is not translationally invariant, it is suitable for 
the study of spherically symmetric boundary-value problems. To demonstrate the utility of 
this approach we compare numerical results obtained from hypercubic and fractal lattices 
with our exact analytical results and show that these results are universal. 

A random walk on a lattice is described by a single probability function, C n t;m , the 
probability that a random walker who begins at location m at t = will be at location n 
at time t. The probability C n t m satisfies a difference equation of generic form: 

C n ,t;m = £ P ^ ~+ n )^-l;m, (1) 

where {cr} n is the set of locations o adjacent to the location n and P(n — > n') represents 
the probability that a walker at the location n will go to the location n' in one step. The 
probability C n t;m obeys the initial condition C n . m = <5 n ,m- Local conservation of probability 
is expressed by S{o-} n -P( n — >• cr) = 1 for all n. 

In this letter we generalize the difference equation ([!]) to allow for creation and anni- 
hilation of random walkers. Specifically, we let random walkers give birth in one location, 
which we label 0, with birth rate a and die in all other locations with uniform death rate z. 
Walkers are created or destroyed in a given location in proportion to the number of walkers 
at that site, where a and z are the constants of proportionality. Of course, a acts as a birth 
rate if a > 1; if a < 1, it is really a death rate. A similar interpretation applies to z. 
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Birth rates and death rates are properties of populations rather than of single individuals. 
Thus, rather than solving an initial-value problem for a single random walker, we study a 
large population of random walkers. We represent such a population by a distribution G n ,t, 
which denotes the number of random walkers at location n at time t. The distribution G n j 
satisfies the same recursion relation as C nji;m except for factors of a and z: 

{ 2E W /(^n)G, ( -i (n^ Wo), 

Gn,t = (2) 

( aP(0 -> n)G ,t-i + z E{^o}» P{° -> n)G CTji _i (n £ Mo). 
Note that G nj t > for all n and t. We are interested in the asymptotic behavior of G n j 

as t — > oo. The specific choice of initial distribution G n0 is unimportant, aside from the 
requirement that it be normalizable, because a random walk is a diffusive (dissipative) 
process and details of G ni0 are irretrievably lost as time evolves; all initial distributions lead 
to the same large-time behavior. We are interested in steady-state solutions of Eq. (Q); the 
existence of such solutions imposes a relationship between the birth and death rates. 

To describe a dissipative system interacting with a hyperspherical boundary we formulate 
our random walk model on a spherically symmetric lattice in D-dimensional space 0. We 
consider the boundary to lie at the origin. (The general case of a hyperspherical boundary of 
arbitrary radius is considered in Ref. JTOf.) We take location n to be the volume bounded by 
two concentric D-dimensional hyperspherical surfaces of radii R n -i and R n . If the random 
walker is at location n at time t, then at time t + 1 the walker must move outward to location 
ri+1 with probability P ut{ n ) or inward to location n — 1 with probability P in (n), where 
Pout{ n ) + -Pin(^) — 1 so that probability is locally conserved. Because a walker at the origin 
(n = 1) must move outward, we take P out (l) = 1 and -Pin(l) = 0. Then, the equation for a 
D- dimensional random walk in Eq. (Q) reduces to a one- dimensional recursion relation: 
' zP hl (n + l)G n+ i jt _i + zP out (n - l)G n _ 1)t _i (n > 3), 

G n>t = | zP in (3)GVi + aGn-M-i ( n = 2 )' ( 3 ) 

k zP m {2)G 2 ^ l {n = 1). 

In Ref. we take the probabilities of moving out or in to be in proportion to the hyper- 
spherical surface areas that can be traversed on each move. Let Sd(R) = 2n D l 2 R D ~ l /T(D /2) 



represent the surface area of a D-dimensional hypersphere. Then, for n > 1, 

P (n) — SoiRn) p / ^ _ Sp^Rn-l) 

Sd(Rti) + Sd(Rti-i) S D (R n ) + SoiRn-l) 

For dimensions other than D — 1 and D = 2, when we take i? n = n, the difference 
equation in (|j) cannot be solved in closed form ||. Thus, we propose that the probabilities 
in Eqs. (|4]) be replaced by bilinear functions of n, which are uniformly good approximations 
to P ut( n ) an d P>n(n) in the range D > when R n = n |J: 

Pont(n) = — — and Pm{n) = . (5) 

v J 2n + D-3 K ' 2n + D-3 K J 

In Ref. @] it was shown that this crucial simplification in Pi n and P ou t preserves the config- 
urational entropy. Now, Eq. (||]) can be solved in closed form for all D > 0: 

r(^)r (n + C -2) , ^ (j) , "^(f,y ; n + ^ ;y ») 

2»-'T(B)r („ + ^2) fc 2niy 1 + (f - l) „F, (1,1;^;^)' 1 

From the solution in Eq. (^|) we can obtain closed-form expressions for spatial and temporal 
moments of the random walk. 

To study critical behavior in the random walk models in Eq. or Eq. ([3]) we take F, 
the fraction of walkers at the boundary (the location where random walkers are born), 
as our order parameter. For the case of polymer adsorption, this observable corresponds to 
the fraction of the polymer adsorbed on the boundary. To construct the fraction F we let 
Nt = J2 n G-n,t be the total number of random walkers at time t. We restrict our attention to 
initial distributions for which iVo is finite so that N t is finite for all t. Then, F t = Go tt /N t 
represents the fraction of all random walkers at location at time t. 

We find that, independent of the choice of G n ,o, the asymptotic behaviors of N t and 
F t are determined by the values of a and z. Specifically, we obtain the generic (lattice- 
independent) result that the positive quadrant of the (a, z) plane is partitioned into four 
distinct regions by two boundary curves as shown in Fig. [j]. One boundary curve, labeled 
B\ in Fig. [TJ, is a straight line passing through the origin. To the left of B±, F t vanishes as 
t — > oo; to the right of B±, F t approaches a positive finite value as t — >• oo. When the lattice 



dimension D < 2 the equation for the boundary line B\ is z = a; as D increases beyond 2 
the boundary line remains straight but the slope of Bi begins to decrease with increasing 
D. The transition at D = 2 is a reflection of Polya's theorem J7|, which states that when 
D > 2 the probability of an individual random walker returning to a location is less than 
unity. 

The second boundary curve shown in Fig. [I] is labeled B2. The first part of B2 is a 
straight-line segment, z = 1, extending from the z axis to the boundary line B\. This line 
segment connects to the second part of B2, which is a curve that approaches z = as a — > 00. 
The equation for this second part depends on D. [For a D = 1 lattice, where the probabilities 
of moving to the left or right are both |, this curve is given by z = 2a/ {a 2 + 1) (a > 1).] 
Above B 2 , N t — > 00 as t — > 00; below 5 2 , iV t — > as t — > 00. On _B 2 the total number of 
walkers approaches a finite number N(a) as i — > 00. On the curved portion of B 2 , N(a) > 0; 
on the straight-line portion of B 2 , N(a) > for D > 2, while iV(a) = for D < 2. This 
transition at D = 2 is yet another manifestation of Polya's theorem. 

As we cross the boundary line B\, the limiting value of the function F t as t — > 00 is 
continuous. We are particularly interested in crossing from one side of i?i to the other along 
the boundary curve B 2 that divides the upper region, where Nt — > 00, and the lower region, 
where Nt — > as t — > 00. The interpretation of lim^oo Nt being finite and nonzero is that 
the distribution G n t approaches a steady-state, where there is a balance between random 
walkers being created at location and annihilated in all other regions. 

We focus on B 2 because it is only on this curve that a steady state is reached as t — > 00. 
Along B 2 , F(a) = lim^oo F t undergoes a second-order phase transition at the critical point 
(a c ,z c = 1), which is situated at the intersection of B\ and B 2 . On B 2 , F(a) = when 
a < a c (even though the limiting value of Nt may be 0), and both N(a) and F(a) are finite 
positive numbers when a > a c . The curved portion of B 2 is the locus of all points in the 
positive quadrant of the (a, z) plane for which both N(a) and F(a) are finite and nonzero. 

Many of the qualitative features of Fig. [I] can be derived directly from an analysis of 
Eq. @. To determine the boundary line B\ we introduce a change of variable: G n> t = z t H n t . 
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The distribution H n j satisfies the recursion relation for a random walk with birth rate a/z at 
location and no births or deaths at any other location. Let n denote the probability that 
a random walker at location will ever return to location 0. Suppose a random walk with 
birth rate a/z and death rate 1 begins at t — 0. Of the i?o,o walkers who begin at location 
0, only a fraction n of them will return to location to give birth to new walkers. Of these 
new walkers, again only Ho of them will return to location to give birth again, and so on. 
Hence, to find the total number of walkers who are ever born we must sum a geometric series 
whose geometric ratio is the quantity allo/z. If allo/z < 1, the geometric series converges 
and the total number of random walkers ever born is finite. As time t increases, the walkers 
diffuse away from location 0. Thus, as t — > oo, the ratio 

p _ (*a,t Ho t 

En G n j J2 n H n ,t 

vanishes. On the other hand, if allo/z > 1, both Ho t and J2 n H n ,t diverge at the same rate, 
and the ratio F t approaches a nonzero limit (that lies between and 1) as t — > oo. 

The transition between F t — > and F t — > finite limit occurs on the line z = alio. This 
is the equation of the boundary line B\. Polya's theorem states that for any random walk 
n = 1 when D < 2 and il < 1 when D > 2. This theorem explains the transition in the 
slope of the line B\ at D = 2 0. In particular, since z c = 1 for all D > 0, we find a c = l/Ho- 

The shape of the curved part of B 2 depends on the choice of lattice; it is not universal. 
However, the straight-line portion of B 2 is universal and is easy to understand. Points (a, z) 
such that a < a c and z near 1 lie to the left of B\. Thus, F t becomes vanishingly small as 
t — > oo. Hence, the effect of the birth rate a on the total number of walkers is negligible. 
The growth or decay of the total number of walkers depends only on z; if z < 1 then iV t — ► 
as t — > oo, and if z > 1 then N t —>■ oo as t — > oo. 

On the straight-line portion of the curve B 2 , where a < a c and z — 1, the limiting value 
of Nt depends on the dimension D. If D < 2 then a c — 1. Thus, on this portion of B 2 a 
fraction 1 — a of random walkers who arrive at location at a given time step must die at 
the next time step. But by Polya's theorem all random walkers visit location repeatedly. 
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Hence, the total number of random walkers N t must vanish as t — > oo. On the other hand, 
if D > 2 we have n < 1. Thus, the fraction 1 — n of random walkers who originate in 
location never return to 0. These random walkers never die because z — 1. Hence, Nt 
approaches a finite positive number as t — > oo. 

There is a change in the form of the transition at D — 4. When D < 4 the slope of 
B 2 is continuous. However, when D > 4 an elbow appears in £?2 at the critical value a c . 
Specifically, when D > 4 the slope of B 2 is for < a < a c ; just above a c the slope abruptly 
becomes — Uo/(T — 1). Here, T is the expected time for a random walker who begins at 
location to return to when a = 1 and z — 1. A quick heuristic argument yields this 
result HJ]: When D < 4 we know that T, the first temporal moment of Co,t;0> is divergent; T 
is finite when D > 4. In a steady state all random walkers at location leave this location 
and in a z = 1 model only the fraction n ever return. The random walkers who return to 
location do so in T steps on average. These returning random walkers experience a death 
rate z for T — 1 of these T steps. Thus, the expected number of random walkers who actually 
return to location is decreased by the factor z T ~ x . Hence, after T steps we expect to find 
the number of walkers at location to be multiplied by the factor aU. z T ^ 1 . The condition 
that there be a steady state is therefore an -2 T_1 = 1. Near the critical point a = a c + 5 and 
z = 1 — e as 5, e — > + . To first order in 5 and e the steady-state condition gives the slope 
of B 2 as a —>■ a+. Note that this argument is valid only for a > a c . 

While the above argument is only valid in the neighborhood of a c , the same reasoning 
yields the entire curve z(a) in the limit D — > 00. In this limit T = 2. Hence, the steady-state 
condition gives z = a c /a. 

To determine the nature of the phase transition in F(a) along B 2 , we let g n = lim^oo G n> t- 
The steady-state distribution g n satisfies the difference equation 



[ aP(0 -> n)g + z E{^o} n P{p -> n)g ff (n e {a} ), 
which is the time-independent version of Eq. (Q). For a steady-state solution having a finite 
number of random walkers the sum J2n 9n exists. If we sum both sides of Eq. ([?]) over all 
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locations, we obtain 



J2dn = (a- z)g + z^g n . (8) 

n n 

Since the sum J2 n 9n is nonzero we obtain a formula for F(a) in terms of a and z(a): 

F(a) = -*>- = ^44 - 0) 

The result in Eq. (||) is valid on the curved part of B2; it is also valid on the straight-line 
portion of B2 when D > 2. On B2 we must treat z as a function of a. We emphasize this 
dependence by writing z(a) and by treating the fraction F as a function of a only. 

We have studied Eq. (|7|) for several types of lattices. We have determined that along the 
curve B 2 the steady-state distribution fraction F(a) behaves like 

F(a) ~ C(D)(a-a c y (a -> a+ D ^ 2, 4). (10) 

The multiplicative constant C(-D) depends on the specific choice of lattice. However, the 
critical exponent v is universal and only depends on the dimension D. Using the probabilities 
in Eqs. (|5|) for the hyperspherical lattice, we derive in Ref. || 

r_o_ (0<J D<2), 
jh (2<D<4), (11) 
1 (D > 4). 

We have verified the universality of these results numerically for more complicated lattices. 
We have computed the critical exponent v for three types of random walks that are not 
analytically tractable: (i) a spherically symmetric random walk using the probabilities in 
Eqs. (|j) for various values of D; (ii) a conventional random walk on a D = 3 cubic lattice; (iii) 
a random walk on a recursively defined fractal lattice (Sierpinski carpet) || with dimension 
D = ^f- The agreement between the predicted value of the critical exponent v in Eq. flllD 
and that obtained by computer simulation is excellent, which provides strong evidence that 
the second-order phase transition at a c is universal. 
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There is no critical exponent for the special cases D = 2,4; for the case of a spherically 
symmetric random walk using the probabilities in Eqs. fl5|) we find that [|J as a — ► a+ 

( constant 2 



F(a) ~ ^ 



P = 2), 
(£> = 4). 



a — a c 



12 log f- 3 ^- 

This work can be applied to the study of polymer growth in the vicinity of an attractive 
boundary ||. For a hyperspherical boundary of radius m, this phenomenon is described by 
a partial difference equation of the same form as that in Eq. (0) [Ej. Let P(k) represent the 
fraction of a polymer that is adsorbed as a function of the attractive boundary potential k. 
This fraction is analogous to F(a) in Eq. The potential k is a monotonic function of 
the birth rate a, and near the critical point 



P{k) oc „ C(£))(o - a c ) v ' '. (13) 



We find that u' — p (D < 2) and z/ = 1/ — 1 (£) > 2), where z/ is given in Eq. ([ill). Thus, the 
polymer adsorption fraction P(k) exhibits a first-order transition for D > 4 and a tricritical 
point with logarithmic scaling at D = 4. This behavior for P(/c) is illustrated in Fig. 0. 

The analysis in this paper allows us to make a physical prediction for systems where 
excluded-volume effects can be neglected: The adsorption fraction for a solution of polymers 
growing near an approximately spherical attractive boundary (such as a cell membrane) 
exhibits a second-order phase transition with linear scaling. 
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of Energy for support under grant number DE-AC02-76-CH00016. 
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FIGURES 

FIG. 1. Generic phase diagram for the (a, z) plane. Shown on the diagram are the boundary 
curves B± and E>2- To the left of B\ and on B\ the fraction of random walkers at location 0, Ft, 
approaches as t — > oo; to the right of B\ this fraction approaches a finite positive number as 
t — > oo. Above B2 the total number of random walkers, N t diverges as t — > 00; below B2 the total 
number of walkers approaches as t — > 00. On £?2 the distribution of random walkers approaches 
a steady state as t — > 00. The critical point (a c , z c = 1) lies at the intersection of B\ and B2. 

FIG. 2. Plot of the adsorption fraction P(k). For increasing D < 2 the scaling exponent 
increases and the transition becomes weaker until at D = 2 exponential scaling is obtained. For 
increasing D > 2 the scaling exponent decreases and the transition becomes stronger again; this 
is compensated by an increase in the critical binding potential k that is required to bring about 
the transition. At D = 4 we observe a tricritical point with logarithmic scaling. For D > 4 
the transition is first order, indicated by a discontinuity (green shaded region) in P(k) across the 
critical point. 
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